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' Abstract 

Dynamics of confinement, chiral symmetry breaking and thermal 
^ | phase transition is considered at large N c . It is argued that these phe- 

nomena are quantitatively well described within the Gaussian stochas- 
00 ■ tic model of the QCD vacuum. Selfconsistent equations are written 

for the field correlators of the model, yielding important connection 
between gluonic correlation length of the vacuum and the string ten- 
sion. Comparison to other approaches and experimental and lattice 



Oh 



Qh. data is given. 



1 Introduction 

X 

Large N c limit first introduced for QCD in [1] has two implications. First, it 
allows to neglect nonplanar perturbative (P) diagrams and reduce nonper- 
tubative (NP) background diagrams to simple expressions; second, it allows 
to establish hierarchy of physical characteristics, which can be compared to 
experiment. E.g. the decay widths of all mesons T n are 0(l/iV c ) while masses 
M n are 0(N®), in experiment on average T n /M n ~ 0.1 which gives an idea 
of the parameter for the real QCD. Also white objects do not interact in 
the leading N c order, which leads to the so-called topological expansion of 
high-energy scattering amplitudes. All dynamical picture of QCD at large 

* Invited Talk at the International Conference on Problems of Quantum field Theory, 
Dubna, Russia, July 1998 



1 



iV c looks selfconsistent and realistic, and there have been many attempts to 
derive it from the first principles [2] . 

Recently a new and universal method was suggested in QCD [3] which al- 
lows to calculate all amplitudes in terms of a set of basic quantities: field cor- 
relators (FC). The simplest approximation uses only the lowest FC, quadratic 
in field strength, and was called the Gaussian Stochastic Model (GSM). Both 
lattice calculations [4] and theoretical estimates [5] show that accuracy of 
GSM is of the same order as of the large N c expansion, and that GSM yields 
the dominant contribution to the most QCD phenomena in the leading large 
N c order. 

The dynamical input of GSM - the quadratic FC - is usually taken in 
a simple exponential form from the lattice data [4]. In the most recent 
paper [6] the large N c limit was used to obtain selfcoupled equations for 
FC, which in principle allow for selfconsistent determination of FC from the 
QCD Lagrangian. Keeping only the lowest FC in these equations one is 
able to connect NP parameters ( string tension, gluon condensate, correla- 
tion length) between themselves and to the P parameter (Aqcd)- 111 this 
way there appears a realistic possibility to have a unified consistent picture 
of QCD dynamics at large N c , described by the only scale parameter as it 
should be. In this talk we discuss confinement (section 2), chiral symme- 
try breaking (CSB) (section 3), thermal phase transition (section 4), with 
concluding remarks in section 5. 

2 Confinement at large N c 

Lattice calculations performed at R < 1.5 fm confirm that static fundamen- 
tal quarks are confined by the linear potential V(R) = aR [7] even in the 
presence of dynamical fermions [8] which in principle should screen the quark 
charges at large distances. Since dynamical quarks break the fundamental 
string in the next 0(l/iV c ) order this means that in the region R < 1.5 fm 
the fundamental string does not break and the leading 0(iV°) approximation 
is valid. 

Another feature of the confinement, also seen in lattice calculations at 
R < 1.5 fm [7] is the so-called Casimir scaling [9], namely that the string 
tension a(j) for static charges in the j representation of the color group 
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SU(iV c ) obey the law 

g(j) C 2 {j) () 

a(fund) C 2 {fundy 1 ' 

where C 2 (j) is the quadratic Casimir operator, 

N 2 - 1 

C 2 (adj,N c )=N c , C 2 (fund,N c ) = -^-. (2) 

The breaking of the adjoint string is the 0(1 /N 2 ) process which happens 
at large distance Rb r (N c ) which grows with N c , as can be seen from the 
expression [9] 

< W adl (C) >= C ie xp(-a adj RT) + ^exp(-V s (R)T). (3) 

c 

Here V S (R) is the screened quark potential. 

We now show that all the above features are nicely described in the GSM. 
To this end one writes the nonabelian Stokes theorem [5,10] 

< W(C) >=< -^-Ptr exp ig [ A^dx^ >= (4) 
<P tr exp ig da^F^u, z ) > 

where we have defined 

rx 

F^(u,z ) = $(z ,w)F^(w)$(w,zo),$(x,?/) = Pexpig / A^dzp (5) 

J y 

and integration in (4) is over the surface S inside the contour C, while z is 
an arbitrary point, on which < W(C) > evidently does not depend. In the 
Abelian case the parallel transporters Q(zo,u) and <!>(u,zo) cancel and one 
obtains the usual Stokes theorem. 

Note that the nonabelian Stokes theorem, eq. (4), is gauge invariant even 
before averaging over all vacuum configurations - the latter is implied by the 
angular brackets in (4). 

One can now use the cluster expansion theorem to express the r.h.s. of 
(4) in terms of FC, namely [3,5] 

< w{c) >= ^r ex pE ^f" / <Mi)<M2)-<M«) < ^(i)-^) > (6) 



where lower indices of da^ v and F^ v are suppressed and F(k) = F^ kUk (u^ k \ z ). 

Note an important simplification - the averages <C F(l)...F(n) ^> in the 
color symmetric vacuum are proportional to the unit matrix in color space, 
and the ordering operator P is not needed any more. 

Eq. (6) expresses Wilson loop in terms of gauge invariant FC, also called 
cumulants , defined in terms of FC as follows: 

< F(1)F(2) >=< F(1)F(2) > - < F(l) >< F{2) > (7) 

< F(1)F(2)F(3) >=< F(1)F(2)F(3) > - < F(1)F(2) X F(3) > - 

- < F(l) >< F(2)F(3) > - < F{2) >< F(1)F(3) > - < F(l) >< F(2) >< F(3) > 

In the lowest approximation, which corresponds to the GSM, one keeps 
only the quadratic in F term, namely 

D^vxo = < ^u(x)<^(x, y)F Xu (y)<&(y, x) > (8) 

The form (8) has a general decomposition in terms of two Lorentz scalar 
functions D{x — y) and D±(x — y) [3] 

D^\a = (S^x ~ &imt&v>)'D{x ~y)+ (9) 
+ \d^{[{hx ■ 5 ua - hJ uX ) + - y)} 

Here the ellipsis implies terms obtained by permutation of indices. It is 
important that the second term on the r.h.s. of (9) is a full derivative by 
construction. 

Insertion of (9) into (6) yields the area law of Wilson loop with the string 
tension er 

< W(C) >= exp(-aS min ) 

a= 1 -jD(x)d*x(l + 0(FT*)) (10) 

where O(FTg) stands for the contribution of higher cumulants, and S min is 
the minimal area for contour C. 

The string tension for heavy (static) quarks is an infinite sum of FC of 
the field F 14 = E 1 integrated over the plane 14: 

n Ti—1 

° ~ E —\ II d ^ < E 1 (0)E 1 (r 1 )E 1 (r 1 + r 2 )...E 1 Q2 r ) » 

n n - i 
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One can identify parameter of expansion in the sum above to be (only even 
powers of n enter the sum) 

C = (E.T 2 ) 2 

where T g is the gluonic correlation length in the vacuum and E\ = g 2 < 
(E®) 2 >= ~ (0.2GeV) 2 where G 2 is the standard gluon condensate. 

Naively one would expect that the correlation length T g is of the order of 
confinement radius R C) T g ~ R c ~ Aq^ d , in which case since E\ ~ ^qcd the 
parameter ( is ( ~ 1 and all FC are equally important (that would be true 
for ordinary FC, but connected FC may have additional small parameter at 
large n due to cancellation in (7)). 

However lattice calculations confirm that T g is much smaller [4], indeed 
T g 0.2 4- 0.3/m and therefore parameter ( is small 

C = 0.04 4-0.1 

The regime ( <C 1, which seems to be characteristic for real QCD, can 
be called the regime of the weak confinement . In this case the dynamics 
of quarks and gluons is adequately described in most cases by the lowest 
(Gaussian) correlator [3,4]. 

Note that D 1 does not enter a, but gives rise to the perimeter term and 
higher order curvature terms. On the other hand the lowest order perturba- 
tive QCD contributes to Di and not to D, namely the one-gluon-exchange 
contribution is 

D?rt{x) = (n) 

Nonperturbative parts of D(x) and D 1 (x) have been computed on the 
lattice [4] using the cooling method, which suppresses perturbative fluctua- 
tions. As one can see both functions are well described by an exponent in the 
measured region, and D 1 (x) ~ \D(x) ~ exp(—x/T g ) , where T g ~ 0.2 fm. 

The string tension (10) can be computed from the lattice data [4] extrap- 
olated to small distances, and agrees within 10-20% with the standard value 
a fa 0.2GeV 2 . Hence the Gaussian correlator alone can explain the string 
tension. 

We conclude this chapter with discussion of confinement for charges in 
higher representations. As it was stated in the previous chapter, our defini- 
tion of confinement based on lattice data, predicts the linear potential be- 
tween static charges in any representation, with string tension proportional 
to the quadratic Casimir operator. 
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Consider therefore the Wilson loop (4) for the charge in some represen- 
tation; the latter was not specified above in all eqs. leading to (10). One can 
write in general 

A^x) = A a ^T a , tr(T a T h ) = \d ab (12) 

Similarly to (6) one has for the representation j = (m 1? m 2 ...) of the group 
SU(N) with dimension N(j) 

1 00 (?a) n r 

< W ( C ) = TT7^ tr J ex P E ^f- / da(l)...da(n) < F(l)...F(n) > (13) 

and by the usual arguments one has Eq.(9) . 

Due to the color neutrality of the vacuum each cumulant is proportional 
to the unit matrix in the color space, e.g. for the lowest cumulant one has 

< F(1)F(2) > ab =< F c (l)F d (2) > T c an T d ah = (14) 
=< F e (l)F e (2) > -J— Tl n T c nh = A^C 2 (j) ■ l ab , 

c 

where we have used the definition 

T c T c = C 2 (j)l (15) 
and introduced a constant not depending on representation, 

A (2) = < F e (l)F e (2) >, (16) 

c 

and also used the color neutrality of the vacuum, 

< F e (l)F e (2) > 

For the next - quartic cumulant one has 

< F(1)F(2)F(3)F(4) > a£ =< F ai (l)F a2 (2)F a3 (3)F a4 (4) > x (18) 

x7$7$7$7£ = A?(C 2 (j)) 2 S a e + A- { 2\T a ^T^T a % £ 

Thus one can see in the quartic cumulant a higher order of quadratic Casimir 
and higher Casimir operators. 



< F c (l)F d (2) >= 5 cd ^_ j '- (17) 
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The string tension for the representation j is the coefficient of the diagonal 
element in (14) and (18) 

<J) = C 2 (j) J 9 -^d*x + 0(Cl(j)) (19) 

where the term 0(Cf (j)) contains higher degrees of C^j) and higher Casimir 
operators. 

Comparing our result (19) with lattice data [7] one can see that the first 
quadratic cumulant should be dominant as it ensures proportionality of a(j) 
to the quadratic Casimir operator. 

Thus we see that GSM yields a simple confinement picture which is con- 
sistent with lattice data and large N c considerations at least in the region 
R < 1.5 fm. At larger R and fixed iV c and for adjoint charges the screening 
occurs which needs higher cumulants, and this happens in the 0(1/ N%) order. 

3 CSB at large N c 

In the large N c limit the phenomenon of CSB is supported by the Coleman- 
Witten theorem [11], whereas at N c = 2, 3 lattice data show clearly CSB in 
the order parameters of quark condensate < ^ (0)^(0) >, which is nonzero 
for T < T c . Also parity doublets are missing in the hadronic spectrum. 
There is still another important feature of CSB seen in the heavy-light qQ 
system, namely the scalar confining interaction for the light quark, which 
clearly signals CSB. In the present section we shall present results of recent 
studies of this system [12,13], which is the simplest from the dynamical point 
of view. 

We start with the quark Green's function and write the effective quark 
Lagrangian in presence of a static source, using the averaging of the partition 
function over gluonic field A^. 

To take into account the static source we consider the generalized coor- 
dinate gauge [10] and express through F^ v as 

where the contour C starts at x^ and is described by z^(x, s) (in the usual 
coordinate gauge z^XjS) = 3x^,0 < s < 1). The effective Lagrangian is (a 
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more extended version of this derivation see in [12]). 




x 



(21) 



where we have defined 



J^'(Z,W)= du a dv^Sa-fSps ~ SasSfa) 




X 



x D(u - v) 

and D{u) is defined in (9) yielding the string tension 



(22) 



a = - / d 2 uD(u) 

2 J— oo 



(23) 



Note that we have neglected in (21) higher field correlators, which were 
argued above to yield subdominant contribution. 

The Lagrangian (21) can be used to obtain equations for the quark 
Green's function S in the large N c limit, where the following rule of re- 
placement holds 



One obtains a system of equations for the quark Green's function S and the 
mass operator M 



The system of equations (25,26) is exact in the large N c limit, when higher 
correlators are neglected and defines unambiguously both the interaction 
kernel M and the Green's function S. One should stress at this point again 
that both S and M are not the one-particle operators but rather two-particle 
operators, with the role of the second particle played by the static source. It 
is due to this property, that S and M are gauge invariant operators, which 
is very important to take confinement into account properly. Had we worked 



(24) 
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with one-particle operators, as is the habit in QED and sometimes also in 
QCD, then we would immediately loose the gauge invariance and the string, 
and hence confinement. 

The CSB can manifest itself in solutions of (25,26) in several ways. One 
is the appearance of the nonzero chiral condensate 

< $(0)V>(0) >= iN c tr S(0, 0) (27) 

This was estimated using the relativistic WKB method in [12] to be 

<^(oM0)>~-iVJ(^,a 3 / 2 ) (28) 

where T g is the gluonic correlation length [3,4]. 

Another manifestation of CSB is the scalar confinement, which is seen 
at large distances r. Indeed one write expansion for the Green's function 
S and M in the inverse powers of the string mass M str = err + m, and the 
time-averaged Green's function S satisfies an equation [13]: 

[-tfd-i(m + a\z\)]S(z,w) = 5^\z-w) (29) 

In (29) the string term or enters as a scalar, which signals CSB. 

Equations (28) and (29) exemplify the connection between CSB and con- 
finement in the large N c limit. The necessary appearance of CSB in this limit 
was proved earlier in [11] but no hint to the possible mechanism was given. 
Here we demonstrate that CSB occurs due to the string formation (which is 
contained in factor a in (28), (29) and the kernel J in (25), (26)), but this 
effect comes to the existence only due to the solution of nonlinear equations 
(25), (26). One should mention that these equations in contrast to those of 
NJL, are nonlocal and therefore do not need cut-off. The CSB solution is 
not obtained by perturbation expansion of (25), (26), but rather is an extra, 
nonperturbative solution existing due to the nonlinearity. A similar situation 
occurs in the NJL model. 

4 Phase transition at large N c 

To describe the phase transition in QCD we shall use the basic idea, suggested 
in [14], that the confining phase is governed by the NP fields, while in the 
deconfining phase one has P fields in the NP background of magnetic fields. 
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Therefore one should first introduce the background formalism where both 
P and NP fields enter. 

We derive here basic formulas for the partition function, free energy and 
Green's function in the NP background formalism at T > [15]. The total 
gluonic field is split into a perturbative part a M and NP background 

A^ = B^ + a„ (30) 

where both and are subject to periodic boundary conditions. The 
principle of this separation is immaterial for our purposes here, and one can 
average over fields B^ and a M independently using the 'tHooft's identity^ 

=<< exp(—S(B + a) > a >B 

with arbitrary weight f]{B). In our case we choose rj(B) to fix field correlators 
and string tension at their observed values. 
The partition function can be written as 

Z(y,T,fi, = 0) =< Z(B) > B , 

Z(B) = N J D<Pexp(- dr J d 3 xL(x,r)) (32) 

where <fi denotes all set of fields a M , N is a normalization constant, 

and the sign <>b means some averaging over (nonperturbative) background 
fields .Bp. 

The inverse gluon propagator in the background gauge is 

Wf v = -D\B) ab ■ 5, v - 2gF^ v {B)r b (33) 

where 

(D x ) ca = dxSca - igT h ca B\ = d x 5 ca - gf bca B b x (34) 
Integration over ghost and gluon degrees of freedom in (32) yields 

Z(B) = N'(detW(B)); e l J 2 [det(-D,(B)D,(B + a))] a=fj x 

oo cl s: I 

x{l + E -f (a = jj)}exp{--JW- l J)j^ d^f^b) (35) 



1 private communication to the author, December 1993. 
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One can consider strong background fields, so that gB^ is large (as com- 
pared to Aq CD ), while a s = in that strong background is small at all 
distances. 

In this case Eq. (35) is a perturbative sum in powers of g n , arising from 
expansion in {ga^) n . 

In what follows we shall discuss the Feynman graphs for the free energy 
F(T), connected to Z[B) via 

F(T) = -Tin < Z{B) > B (36) 

We are now in position to make expansion of Z and F in powers of ga^ (i.e. 
perturbative expansion in a s ), and the leading-nonperturbative term Z ,F 
- can be represented as a sum of contributions with different N c behaviour of 
which we systematically will keep the leading terms 0(N 2 ), 0(N C ) and 0(N®). 

To describe the temperature phase transition one should specify phases 
and compute free energy. For the confining phase to lowest order in a s 
free energy is given by Eq.(36) plus contribution of energy density e at zero 
temperature 

7r 2 V (2m T) 3 / 2 

F(l) = sV 3 - -V 3 T 4 -T]P e-"Wr + 0(1/7V C ) (37) 

where e is defined by scale anomaly [16] 

1 1 n 

£ ~ ~Y Nc 32^ < ^ B))2 > (38) 
and the next terms in (37) correspond to the contribution of mesons (we keep 
only pion gas) and glueballs. Note that e = 0(N 2 ) while two other terms in 
(37) are 0(iV c ). 

For the second phase (to be the high temperature phase) we make an 
assumption that there all color magnetic field correlators are the same as in 
the first phase, while all color electric fields vanish. Since at T = color- 
magnetic correlators (CMC) and color-electric correlators (CEC) are equal 
due to the Euclidean 0(4) invariance, one has 

< {^ U (B)) 2 > = < (F^ u ) 2 > e i + < (F^ u ) 2 > magn ; < F 2 > ma gn = < F 2 > e ; 

(39) 

The string tension a which characterizes confinement is due to the electric 
fields [5], e.g. in the plane (i4) 

2 

cr = o-e = y J j d 2 x< trE i (x)$(x,0)E i (0)(f)(0,x) > +... (40) 
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where dots imply higher order terms in E^. 

Vanishing of <je liberates gluons and quarks, which will contribute to the 
free energy in the deconfmed phase their closed loop terms with all possible 
windings. As a result one has for the high-temperature phase (phase 2). 

F(2) = \eV, - {Nl - 1)V S ^- - ^N C V 3 T% + 0(N° C ) (41) 

Comparing (37) and (41), F(l) = F{2) at T = T c , one finds in the order 
0(iV c ), disregarding all meson and glueball contributions 

Tc = l¥W-i) + ¥^»/J (42) 

For standard value of G 2 = ^ < F 2 >= O.OUGeV 4 (note that for n f = 
one should use approximately 3 times larger value of G 2 ) one has for SU (3) 
and different values of n f = 0,2,4 respectively T c = 240,150,134 MeV. 
This should be compared with lattice data [16] T c (lattice) = 240, 146, 131 
MeV. Agreement is quite good. Note that at large N c one has T c = 0(iV ( ?) 
i.e. the resulting value of T c doesn't depend on N c in this limit. Hadron 
contributions to T c are 0(iV~ 2 ) and therefore suppressed if T c is below the 
Hagedorn temperature as it typically happens in string theory estimates. 



5 Conclusion 

In all examples considered above gluon correlators and condensates entered 
as a given input. It is important to find equations for FC which define 
them up to an overall scale to Aq C d- These equations have been suggested 
in [6] and are derived on the same physical basis as for the quark Green's 
function in section 3. Namely one assumes that a gluon is propagating in the 
nonperturbative background, described by FC, and one obtains equations 
for the gluon Green's function (in the field of the static charge, so that the 
Green's function is gauge-invariant, as in the case of the quark in section 3). 

J M^ p \x,z)G pu (z,y)d 4 z = 5^(x-y), (43) 

where the mass operator is approximately equal at large distances to 
M (2,2) , where 

M$ 2 \x,y) = ^ 4 \x - y)[J xx (x,y)5, u - J»„(x,y)\, (44) 
G 2 
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For more extended treatment and derivation of equations (47), (48) the reader 
is referred to [6]. 

The kernel Ja m in (44) is expressed through the field correlator 
< FF >. To make equations selfconsistent one should express the latter 
through the gluon Green's function G. This is possible since one can always 
refer the color indices in each term of J2 a =i < F a F a > to the group of fields 
b^. Then one can write symbolically 



where ellipsis stands for higher cumulants. 

The system of equations (43-44) allows for a nonperturbative solution, 
which violates the scale invariance present in the equations. This solution 
is defined by fixing one nonperturbative scale, e.g. the string tension a. 
Then equations (43-44) predict that i) both field correlators D(x), Di(x) [3] 
exponentially decay at large x; D(x), Di(x) ~ exp(—x/T g ) in agreement with 
lattice data [4], and ii) the gluon correlation length T g is connected to a as 



Insertion of the standard value a « 0.2GeV 2 yields T g w 0.2/m, which is in 
good agreement with lattice data [4] . 

We have derived equations for the gluon and quark propagators in the field 
of a static source. These equations possess symmetry (chiral for the quark 
and scale invariance for the gluon) which is violated by the nonperturbative 
solutions. One obtains in this way the CSB due to the confining kernel, and 
the confining kernel itself satisfies nonlinear equations. Properties of this 
kernel are in agreement with lattice data. 

Thus the Gaussian model using the simplest FC can describe both qual- 
itatively and quantitatively the basic QCD phenomena at large N c . 



This work was supported by the RFFI grants 96-02-19184a, 96-1596740 
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